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ABSTRACT 

There have been several recent claims of black hole binaries in globular clus- 
ters. I show that these candidate systems could instead be ultracompact X-ray 
binaries (UCXBs) in which a neutron star accretes from a white dwarf. They 
would represent a slightly earlier evolutionary stage of known globular cluster 
UCXBs such as 4U 1820-30, with white dwarf masses ~ O.2M , and orbital pe- 
riods below 5 minutes. Accretion is slightly super-Eddington, and makes these 
systems ultraluminous sources (ULXs) with rather mild beaming factors b ~ 0.3. 
Their theoretical luminosity function flattens slightly just above Z/Edd and then 
steepens at ~ 3LEdd- It predicts of order 2 detections in elliptical galaxies such 
as NGC 4472, as observed. 

The very bright X-ray source HLX-1 lies off the plane of its host SOa galaxy. 
If this is an indication of globular cluster membership, it could conceivably be a 
more extreme example of a UCXB with white dwarf mass M2 — 0.34M Q . The 
beaming here is tighter (b ~ 2.5 — 9 x 10~ 3 ), but the system's distance of 95 
Mpc easily eliminates any need to invoke improbable alignment of the beam for 
detection. If its position instead indicates membership of a satellite dwarf galaxy, 
HLX-1 could have a much higher accretor mass ~ 1000M Q 

Subject headings: accretion, accretion disks - X-rays: binaries - X-rays: galaxies 
- globular clusters 



Introduction 



Globular clusters (GCs) host 10 - 20% of the X-ray sources in galaxies like our own. 
This fraction rises to < 50% in early-type galaxies (e.g. Sarazin et al., 2003). In all galaxies 
this implies a far higher incidence of X-ray sources per unit mass in GCs than in the field - 
in the Milky Way by about a factor 100 - 1000. The reason for this is of course that GCs can 
form close binaries by dynamical capture, and many of these evolve into low-mass X-ray 
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binaries (LMXBs). In early-type galaxies this process continues long after star formation in 
the field has declined, making the GC sources more prominent. 

A sign of the unusual formation channels available in GCs is the prominence there of 
ultracompact X-ray binaries (UCXBs) - systems with very short orbital periods P down to 
~ 10 minutes. For example, 5 out of the 7 LMXBs with P < 30 minutes in the 2011 update 
to the catalogue of Ritter & Kolb (2003) are known to be members of globular clusters. 
These systems all appear to involve a neutron star accreting from a white dwarf companion 
star. In contrast, systems containing black holes are rare or absent from GCs. Clearly 
one would not expect high-mass X-ray binaries similar to Cygnus X-l, but there is so far 
no clear detection in a GC of the soft X-ray transients which constitute the majority of 
Galactic black-hole systems. It is often suggested (e.g. Kalogera, King & Rasio, 2004) that 
most black holes are dynamically expelled from GCs because the mass contrast with the 
other stellar populations makes them vulnerable to the Spitzer mass-segregation instability 
(Spitzer, 1969). 

There have nevertheless been a number of claimed or suggested detections of black-hole 
X-ray binaries in GCs (Maccarone et al., 2007; Brassington et al., 2010; Shih et al., 2010; 
Maccarone et al., 2011) on the basis of luminosities up to 4.5 x 10 39 erg s" 1 , exceeding the 
Eddington limit ^Edd — 10 38 erg s" 1 for a neutron star, together with > 2x variability (to 
rule out superpositions of fainter sources). The same criteria applied to the much brighter 
(~ 10 41 — 10 42 erg s~ x ) source HLX-1 have led to the suggestion of a black hole mass 

> 5OOM (Farrell et al., 2009; Wiersema, 2010). HLX-1 is not known to be a member of 
a GC: it is associated with an edge-on SOa spiral galaxy ESO 243-39, lying off the galaxy 
plane in the outskirts of its bulge, and so could be a member of a GC or other type of star 
cluster, or of a satellite dwarf galaxy. 

In the absence of clean dynamical masses (> 3M Q for the claimed GC black holes, and 

> 500M Q for HLX-1) the nature of these sources is not definitively settled. It is reasonable 
to ask if other interpretations are possible, and I attempt this here. 

2. Bright X-ray Sources in Globular Clusters 

The defining feature of the GC black hole candidates is their luminosity, which is super- 
Eddington for a neutron star. Any model of them must therefore allow for high mass 
transfer rates. Since all the potential donor stars in GC binaries have low masses < 1M Q 
this obviously also explains the short lifetimes of the black hole candidates and thus their 
rarity. We have already noted that the types of black hole X-ray binaries found in the field 
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- high-mass systems, and soft X-ray transients, are not found in GCs. This leaves one fairly 
obvious alternative, the ultracompact systems mentioned above. These systems evolve in a 
very simple way, which naturally implies a very high mass transfer rate at a certain epoch 
(see e.g. King, 1988 for a discussion). The condition that a white dwarf donor should fill its 
Roche lobe leads to a relation 

P ~ O.Om^ 1 minutes (1) 

between donor mass M 2 = m 2 M Q and orbital period P (with the white dwarf taken as 
a hydrogen-depleted n = 3/2 polytrope with radius-mass relation R 2 oc M 2 ). These 
very short orbital periods imply that mass transfer in these UCXB systems is driven by 
gravitational radiation, which gives a mass transfer rate 

- M 2 ~ 1 x Kr 3 mJ /3 mJ 4/3 Moyr" 1 (2) 

where M\ = miM Q is the accretor mass. The brightest low-mass X-ray binary in a Galactic 
globular cluster is 4U1820-30 in NGC 6624, with L ~ (4 - 7) x 10 37 erg s -1 . This has an orbital 
period P = 11.4 minutes (Stella, Priedhorsky & White, 1987; van der Klis et al., 1993), and 
so from ([TJ), presumably a current donor mass M 2 — O.O8M . Equation (T5]) then gives 
— M 2 ~ 7 x 10~ 9 M Q yr _1 , in good agreement with the observed luminosity if the accretor 
is a neutron star of ~ 1.4M . The observed X-ray bursts from this system confirm this 
identification, and suggest that the donor is a He white dwarf (Bildsten, 1995; Strohmayer 
& Brown, 2002; Cumming, 2003). I note that Bildsten & Deloye (2004) consider the effects 
of chemical composition and finite entropy and derive slightly more complex relations than 
([U ED; in particular an exponent closer to 4.13 for m 2 in ([21 . However their Fig.l shows 
that the deviations are small enough that the poloytropic approximation is adequate for the 
purposes of this paper, particularly for larger values of the white dwarf mass than considered 
by Bildsten & Deloye. 

The important result here is that the mass transfer rate could have been much higher 
in the past, when M 2 was larger. The initial white dwarf mass in a UCXB is set by its 
prior evolution. UCXBs presumably result from some kind of common envelope evolution 
following the dynamical capture by the neutron star primary of an evolved companion star. 
The degenerate core mass of the companion can range from ~ 0.1 — O.4M for helium, 
and up to ~ O.6M for carbon/oxygen, depending sensitively on how far the companion has 
evolved at the epoch of capture. It is clear that mass transfer rates well above the Eddington 
value are easily possible in UCXBs. Note that LecM * s twice the usual value since we expect 
hydrogen-poor accretion here, i.e. 

L Edd ~ 3.5 x 10 38 (m!/1.4)) ergs" 1 . (3) 

This implies an Eddington accretion rate 

MEdd-SxlO^Meyr- 1 . (4) 
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3. Ultraluminous X-ray Sources in Globular Clusters 

Compact objects accreting above their Eddington rates probably appear as ultralumi- 
nous X-ray sources (ULXs), and it is now generally accepted that a large fraction of ULXs 
are of this type. Disc accretion in this regime was first described by Shakura & Sunyaev 
(1973). In their picture, radiation pressure becomes important at the spherization radius 
R sph ~ 27mi? s /4, where rh is the local accretion rate in units of the Eddington value, and 
R s = 2GMi/c 2 is the Schwarzschild radius of the accretor. Shakura & Sunyaev explicitly 
considered only black hole accretors, but their picture also applies to other accretors provided 
that rh is sufficiently large to ensure i? sp h > accretor radius (this always holds for neutron 
stars with rh > 1 for example). Inside R sp h the disc remains close to the local radiation 
pressure limit and blows gas away so that the accretion rate decreases with disc radius as 
M(R) ~ M(i?/i? sp h) — M^ddiR/Rs)- The disc wind has the local escape velocity at each 
radius, so mass conservation shows that the wind is dense near i? sp h and tenuous near the 
inner disc edge. The centrifugal barrier along the disc axis creates a vacuum funnel through 
which the luminosity escapes. 

Thus super-Eddington accretion produces large apparent X-ray luminosities through 
two effects of super-Eddington accretion (Begelman et al., 2006; Poutanen et al., 2007). 
First, the bolometric luminosity is larger than the usual Eddington limit by a factor ~ 
(1 + lnm). Second, the luminosity of a ULX is collimated by a beaming factor b via scattering 
off the walls of the central funnel, so that the apparent luminosity (inferred by assuming 
isotropy) is 

L sph ^^(l + lnm). (5) 

Eddington ratios rh » 1 producing ULXs can arise in a number of ways in close binary 
accretion. Thermal-timescale mass transfer from a massive radiative donor to a lower-mass 
compact object can give rh ~ 5000 or more, as in SS433 (King, Taam & Begelman, 2000; 
Begelman, King & Pringle, 2006), and nuclear expansion can give similar rates (Rappaport, 
Podsiadlowski & Pfahl, 2005). Super-Eddington rates can arise during disc instabilities 
(King, 2002), and I suggest here that UCXBs offer another steady channel. By contrast, a 
similar Eddington ratio rh » 1 is very difficult to produce in accretion in active galactic 
nuclei. If the AGN sits in a spheroid of velocity dispersion o and gas fraction f g , even 
the extreme assumption of a full dynamical accretion rate ~ f g cr 3 /G, given by suddenly 
removing centrifugal support from orbiting gas, produces only rh < 40 (King, 2010, eqn 6). 
Binaries do better than AGN because the self-gravity of the donor star allows a large gas 
mass to get close to the accretor simultaneously by spiralling in through near-circular orbits. 
(Dynamical disruption of stars in AGN does not achieve this, as the likely parabolic orbit 
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means that the arrival time of the gas is spread out by factors (MsMBH/^star) 1 ^ 2 > 10 3 — 10 4 , 
cf Lodato et a!., 2009.) Accordingly there are no AGN analogues of ULXs. 

Recently, King (2009) noted that the observed inverse correlation between soft X-ray 
emission and blackbody temperature (L soit oc T^ 3 ' 5 ) for some ULXs (Kajava & Poutanen, 
2009) could be understood if the beaming factor varies as 

73 

, (6) 

a form which also follows from simple geometrical arguments about the funnel opening 
angle. (We note that this suggests that significant beaming occurs only for sufficiently 
super-Eddington accretors rh > 8.5.) Mainieri et al (2010) used this form of beaming to 
study the luminosity function of a sample of ULXs out to redshifts z = 0.3, finding excellent 
agreement. 

Accordingly, using ([3j |5]) I adopt 

L sph = 4.4 x 10 36 mim 2 (l + lnm) ergs -1 . (7) 

for the apparent UCXB luminosity for rh > 8.5, setting 6=lin(j5])forl<m<8.5. 

A super-Eddington UCXB system with a 1.4M Q neutron star accretor and an Eddington 
factor rh — 15 has a L sp h ~ 5 x 10 39 ergs -1 , large enough to explain all the putative black- 
hole systems in GCs. The beaming factor here is 6 — 0.32. To see if this is reasonable we 
need to work out the luminosity function of these systems. 



4. Luminosity Function 

Bildsten & Deloye (2004) note that the evolution of UCXBs is so rapid that one can 
derive their GC luminosity function (LF) by assuming staeady state conditions. This gives 

Extending the LF to luminosities where beaming is important is straightforward, follow- 
ing Bildsten & Deloye (2004). They note that UCXBs evolve through the observable range 
of luminosities on a timescale far shorter than the age of GCs or any reasonable estimate of 
the time for their birthrate to change. Hence if we know the time t that a UCXB spends 
above a given luminosity, then the cumulative LF N(> L) is directly proportional to t. 

To find t as a function of L I integrate the evolution equation (j2J) assuming << M\, 
as t oc M\~^ ', where I have generalized the exponent of rri2 in from the polytropic value 
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14/3 to a general /3 to enable comparison with Bildsten & Deloye, and also assumed that 
the current donor mass is much smaller than its initial value. This gives t oc (— M 2 )~ < - /3 ~ 1 ^' 3 . 
If there were no beaming, we would now replace — M 2 by L and get iV(> L) oc t oc L~^~^^ 
and so diV/dL oc L 1 ^' 2 for rh < 1 (Bildsten & Deloye, 2001). 

However beaming changes both the connection between — M 2 and L, and also the num- 
ber of systems we can see: we now have L oc b~ x oc rh 2 (neglecting the logarithmic dependence 
on rh, and iV(> L) oc t ■ b oc m~^~ 1 )/ /3 ~ 2 . Together these give N(> L) oc L - ^ -1 ^ 2 ^ -1 , and 

^ oc L 1 / 2 ^ 5 / 2 (9) 
dL 

for rh > 8.5. If indeed there is essentially no beaming for 1 < rh < 8.5 as suggested just 
below eqn ([7]), the logarithmic relation between L and rh implies that are more UCXBs than 
naively expected between L = L Edd and L ~ 3L Edd , so the luminosity function flattens there, 
i.e. 

^ocexp 0.77(l (10) 
dL V L Edd J 

Enforcing continuity at the breaks L = L^dd, 3.14LEdd, the UCXB luminosity function is 

r L77 , exp[0.77(l-/)], 2.93/- 2 - 38 (11) 
over these ranges, where we have taken /3 — 4.3 and set I = L/L E dd- 



5. Candidate GC Black Holes 

Maccarone et al. (2007, 2011) discuss two black hole candidate sources in the elliptical 
galaxy NGC 4472. They are XMMU 122939.7+075333 (L < 4.5 x 10 39 erg s" 1 ) and CXOU 
1229410+075744 (L < 2 x 10 39 erg s" 1 ). The LF of Kim & Fabbiano (2004) for NGC 4472 
shows that there are some ~ 300 X-ray sources with 4 x 10 37 ergs -1 < L < Lb (their Fig.2: 
note that the displayed LF is scaled upwards by a factor 29.51 for display purposes). If 
a substantial fraction of these are UCXBs, as expected, the theoretical LF (OH [TO]) shows 
that the presence of two ultracompact neutron-star ULXs with the quoted luminosities is 
perfectly reasonable. For the more extreme (XMMU) source the white dwarf companion has 
mass — 0.2M©, the orbital period is P — 4.5 minutes, and the beaming factor is b — 0.32. 



6. HLX-1 



Soria et al. (2010) find a red optical counterpart to HLX-1. At the distance to ESO 
243-49 this has a luminosity compatible with a massive globular cluster. On the other 
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hand this luminosity is also compatible with the nucleus of a stripped dwarf galaxy, and in 
addition Wiersema et al., (2010) find a bright Ha line associated with the source. If this 
comes directly from the accretion flow it rules out a UCXB model, although one has still 
to rule out an origin in a nebula around the source - this appears unlikely in a gas-poor 
environment like a globular cluster. In view of these uncertainties I consider a number of 
possibilities. 

If HLX-1 is a UCXB in a globular cluster we can use (j7j) with mi = 1.4 to find an 
estimated Eddington ratio m 170, with a beaming factor b ~ 2.5 x 10 -3 , which would 
require a white dwarf donor mass M2 = O.34M , and an orbital period P ~ 2.5 minutes. 
As a check we can estimate the distance D m[n to the nearest object plausibly observable 
with this amount of beaming (i.e. without requiring an improbable 'aiming' of the beam at 
Earth). This is given by eqn (12) or (15) of King (2009) (but note that the coefficient in 
the final form of eqn. (15) should be 258 and not 660). This gives D mm ~ llN^ 1 ^ Mpc, 
where N (presumably ~ 1) is the number of objects of this type in the host galaxy. This is 
comfortably smaller than the known distance of 95 Mpc, so it is not unreasonable to be in 
the beam of this system. 

If instead of being a globular cluster UCXB, HLX-1 involves 10M Q black hole accreting 
hydrogen-rich material we halve the coefficient on the rhs of ([7]) and set mi = 10. This gives 
an Eddington ratio m = 90, which is quite possible for either a high-mass X-ray binary or 
a soft X-ray transient. The beaming factor is b ~ 9 x 10~ 3 , and the minimum reasonable 
distance is D min ~ lliV -1 / 3 Mpc. Finally, if HLX-1 is actually in a dwarf galaxy interacting 
with ESO 243-39, it is quite possible that the accretor is the central massive black hole of 
this galaxy with a mass > 1OOOM . In this case the interaction with the larger galaxy can 
trigger accretion (King & Dehnen, 2005). Lasota et al. (2011) consider other possibilities of 
this type. 

7. Conclusions 

This paper has shown that the claimed black hole binaries in globular clusters are also 
understandable as ultracompact X-ray binaries in which a neutron star accretes from a 
white dwarf. The high luminosity arises from the higher mass transfer rates expected from 
progenitors of known UCXBs such as 4U 1820-30, with white dwarf masses ~ 0.2M Q . The 
theoretical luminosity function of these systems agrees with the detected number. UCXBs 
have orbital periods of only a few minutes. Detection of such periods would strongly support 
a UCXB model, while measurement of secure dynamical masses would break the ambiguity 
between neutron star and black hole accretors. However if the systems are indeed UCXBs 
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both of these will prove difficult, as the mild beaming induced by super-Eddington accretion 
suggests that the systems are likely to be face-on. Longer superorbital periods may be seen as 
a result of precessions in either case. Similarly, it will be difficult to decide spectroscopically 
whether the accreting material is hydrogen-depleted or not, as one needs to find a line feature 
securely associated with the binary. 

The very bright X-ray source HLX-1 could conceivably be a more extreme example of 
a UCXB (white dwarf mass M 2 ~ O.34M ) if it is indeed a member of a globular cluster. 
However if it is a member of a dwarf satellite galaxy its mass could be ~ 1OOOM . There is 
currently no easy way to eliminate any of these possibilities. 

Research in theoretical astrophysics at Leicester is supported by an STFC Rolling Grant. 
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